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Abstract 

We have developed a theoretical model to analyze the anomalous cooling of radiative quantum 
dot (QD) excitons by THz radiation reported by Yusa et al [Proc. 24th ICPS, 1083 (1998)]. We 
have made three-dimensional (3D) modeling of the strain and the piezoelectric field and calculated 
the 3D density of states of strain induced quantum dots. On the basis of this analysis we have 
developed a spin dependent Monte Carlo model, which describes the carrier dynamics in QD's 
when the intraband relaxation is modulated by THz radiation. We show that THz radiation causes 
resonance transfer of holes from dark to radiative states in strain-induced QD's. The transition 
includes a spatial transfer of holes from the piezoelectric potential mimima to the deformation 
potential minimum. This phenomenon strongly enhances the QD ground state luminescence at the 
expense of the luminescence from higher states. Our model also reproduces the delayed flash of 
QD ground state luminescence, activated by THz radiation even ~ 1 s after the carrier generation. 
Our simulations suggest a more general possibility to cool the radiative exciton subsystem in 
optoelectronic devices. 
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I. INTRODUCTION 



Semiconductor quantum dots (QD's) are commonly called artificial atoms, due to their 

nnn 

man-made single-atom-like spectrum, consisting of well separated spectral lines. [111 . 1151 . 118I | 
Nevertheless, time-resolved spectroscopy of intraband relaxation sequences has proved to be 
much more difficult with QD's than with free atoms. It is, in distinction to free atoms, not 
possible to observe QD intraband relaxation as sequences of well resolved intermediate steps, 
because the intermediate states are very close in energy and involve many simultaneously 
excited electrons. 

Recently Yusa et al used visible light and near infrared (NIR) radiation in combination 
with terahertz (THz) radiation to study intraband carrier relaxation in strain-induced quan- 
tum dots (SIQD's).^] The experiment of Yusa et al, analyzed theoretically in this work, 
has remained to our knowledge the only measurement in which the details of the QD in- 
traband relaxation processes have been successfully studied by modulating the dynamics of 
the carriers, with a radiation field, which is in approximative resonance with the intraband 
transitions. This has prompted us to make a careful analysis of this experiment in order to 
emphasize the potential of this experimental approach and to support a detailed planning 
of further experiments of this kind. 

In this work we analyze the experiments of Yusa et al using a spin-resolved Monte Carlo 
(MC) carrier dynamics model, which is based on new full three-dimensional (3D) strain and 
electronic structure calculations and furthermore on information extracted from previous 
time-resolved and continuous wave (CW)[12|] luminescence measurements, including 
their computational analysis. We show that deep minima of the piezoelectric potential 
(PEP) play a key role in the observed anomalous THz radiation enhancement of the ground 
state luminescence. 

The electronic structure was calculated using continuum elasticity and the multiband 
k ■ p method, which are well suited for the modeling of relatively large SIQD's (the ground 
state wavefunctions have an effective diameter of approximately 15 — 20 nm). 28|] We have 
made a detailed comparison of the theoretical CW luminescence with experiments and pre- 
vious axially symmetric calculations. On the basis of this comparison we build a simplified 
electron and hole structure. This was necessary to make the MC dynamics simulations 
computationally feasible. 
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Our semiempirical carrier dynamics model is based on an intuitive description of the key 
physical processes governing the carrier dynamics of a SIQD. The separate relaxation pro- 
cesses and their governing rate parameters are taken from several previous experiments and 
from computational analyses of these experiments. Our MC model is a generalization 
of previous master equation models [sl, Q, by including the spin of the carriers and by 
accounting for the piezoelectric potential minima. 

We emphasize that fully ab initio calculations of the studied electron-hole system, in- 
cluding for instance correlation effects in a nonperturbative way, is at present not com- 
putationally feasible. However, we hope that our work will prompt for detailed studies of 
correlation effects as well as for advanced studies of the electronic structure using microscopic 
theories. |27t| 



II. ELECTROELASTICITY AND BAND STRUCTURE 



The understanding of SIOD's has until now relied on a simplified axisymmetric description 



of the electronic structure. [2^ In this work we have carried out extensive three-dimensional 
modeling of the electroelastic, electronic and photonic properties of SIQD's. We emphasize 
the role of the PEP, not present in previous two-dimensional models of the carrier dynam- 
ics. Figure [H^b) shows our SIQD model geometry, which is based on transmission electron 
micrographs of SIQD's. [6|] 



Atomic structure 



The elastic strain and the related piezoelectric potential were calculated using the finite 
element method within the continuum elasticity (CE) approximation. This was motivated 
by the large geometry, spanning over ~ 10^ atoms and the possibility to include piezoelectric 
coupling in the model. In the coupled electroelastic field formulas, the vectors of the stress 
T and the electric flux D are related to the strain S and the electric field E as 



23| 



T 




C e 




S 


D 




—e 




-E 



(1) 



The symmetry of the crystal enters Eq. [T] through the elasticity matrix C and the piezo- 
electric matrix e. The dielectric properties of the materials are described by the dielectric 
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matrix e = el^^s- 

The lattice mismatch was introduced into Eq. [T] by expanding the different materials 
isotropically according to their lattice mismatch with respect to the GaAs substrate Aaj = 
{ttj — aGaAs)/cij, where aj is the lattice constants of either InP or InGaAs and acaAs is the 
lattice constant of GaAs. Equation [T] was solved iteratively, [l| using tetrahedral second order 
elements and the CE parameters given in Table [B 

Figure [2] shows the coraputed strain-induced electron and hole potentials in the middle of 
the quantum well (QW).[2J] The total potential is the sum of the band-edge discontinuity, 
deformation and piezoelectric potentials. The elastic strain gives rise to the main deforma- 
tion potential (DP) minima in the center. The PEP gives, on the contrary, rise to large side 
minima and side barriers of holes, located near the edges of the InP island. The effect of 
the PEP on the center minima is, however, small. 



B. Electronic structure 



he electronic structure of the QD's was calculated using the eight-band k-p method. |10l. 
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2l[ | accounting for band edge discontinuity, strain-induced deformation potentials and 
the piezoelectric potential. The eight-band k ■ p Hamiltonian can be divided in a strain-free 
part and a strain part. The strain-free part is parametrized by the band gap Eg, the spin- 
orbit splitting energy A, the optical matrix element Ep, the valence band offset VBO, the 
modified Luttinger parameters 71, 72 and 73, the conduction band effective mass rrtg and 
Kane's parameter B. The strain part, on the other hand, is parametrized by the deformation 
potentials Oc, a„, and dy. Moreover, it contains the six components of the symmetrized 



strain tensor Sxx, £yy, £zz, £xy, £yz, £xz, and the piezoelectric potential ^pje^-jsl 

The eight-band k ■ p Hamiltonian was discretized using the finite difference technique 

□ 

and diagonalized using the implicitly restarted Lanczos method [201]. For more details on the 



electronic structure calculations see Ref. 
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The material parameters used in the eight- 



band k ■ p calculations are listed in Table [TTl The value of the conduction-valence band 
coupling parameter Ep, given in Table [TTl has been reduced in order to avoid spurious states 



fsee Ref. 
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% and references therein). However, the effect of the reduction of Ep on the 
effective masses was negligible (this was verified by testing the model for the pertinent bulk 
semiconductors) . 
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The depth of the PEP minima (+65 meV, see Fig. [2]) of holes exceeds that of the DP 
minimum (+24 meV) by tens of meV. The PEP minima have consequently many confined 
hole states and it is, therefore, energetically more favorable for the holes to populate the 
deepest PEP states before occupying the DP states. The electrons are, however, mainly 
confined in the DP minimum (-108 meV), which over-rules the shallow PEP minima of 
electrons (-15 meV). Spatially separated electron (located on the DP ground state ei) and 
hole (on the PEP level pi) pairs have several decades smaller probability of recombining 
radiatively than an exciton confined to the DP levels ei and hi. 

Figure [3] shows the probability densities of the lowest optically active electron and hole 
states, confined in the DP minima in the center of the SIQD. Figure H] shows the density of 
states (DOS) of all (bright and dark) confined states, including also the hole states confined 
in the PEP minima (dark). The eigenenergies of Fig. [Hare given with respect to the strain- 
free band edge energy of the QW material. 



C. Photoluminescence 



The PL intensity was calculated by summing over the transition rates from all electron 
states to all hole states. The differential photoluminescence rate (rate of spontaneous emis- 
sion of photons) of light with the polarization ex. is given by 

UE) oc 5^ ^ I (vj/,, |«* . p| ME,) U{Ef) P{[E-{E,- Ef)l V} , (2) 

where g{E) is the density of photon states [g{E) oc E'^ for a 3D system, i.e., in the absence 
of a cavity], Ei and Ef are the eigenenergies of the initial electron states I^E'ei) and the final 
hole states \^hf)- The inhomogeneous linewidth broadening of the transition energies Ef^ = 
Ei — E'j, due to the size distribution of the SIQD, was approximated using the Gaussian 
broadening function P(i?, r),[8, [l3] where F is the inhomogeneous linewidth. The Lorentz 
shaped homogenous linewidth broadening (due to the finite life times of the QD eigenstates) 
was neglected, because it has a vanishing effect on the PL, which is also inhomogeneously 
broadened. In the calculations of Fig. [5] we have assumed complete state filling, i.e., the 
Fermi functions of electrons and holes were fe{Ei) = 1 and fh{Ef) = 1, respectively. On 
the basis of the experimentally observed size distribution of InP stressor islands [l^ we have 
estimated F = 6 meV. 
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Figure shows the simulated PL spectrum, averaged over the polarization, for an ensem- 
ble of SIQD's together with an experimental luminescence spectrum. Figure E] shows also the 
effect of the exact InP island geometry and of the PEP on the PL spectrum. All results are 
based on fully three-dimensional simulations. However, in the axisymmetric model we used 
an InP island having the shape of a truncated cone with the same height, the same InP vol- 
ume and the same bottom area as the angular island, shown in Fig. [T](b). We note that the 
PEP affects mainly the PL peaks corresponding to transitions between excited states (QD3, 
QD4, QD5), whereas the angular island geometry (the deviation from the axisymmetry) 
increases the peak-peak energy separation and sharpens the peak profiles. 

The calculation of radiative transition energies and line strengths shows good agree- 
ment with the experimental intensity distribution and energy spacing of the luminescence 
peaks. {12! Each peak consists of many different transitions which after the convolution of 
inhomogeneous broadening form a smooth and continuous spectrum, well reproducing the 
experimental PL. A characterization of the PL peaks, according to the symmetry of the 
involved eigenstates, is difficult due to the weak confinement and small energy spacing 
(1 — 4 meV) of hole states. Only the two lowest peaks (QDl and QD2) are due to recom- 
bination between clearly confined electron and hole states. The higher PL peaks are due to 
recombination between a few confined electron states and several confined or quasiconfined 
hole states (i.e. the hole states are shared between the DP and PEP minima, with their 
energies merging with the energy continuum of the QW). 

III. DYNAMICAL MODEL 

The presence of the PEP minima makes the relaxation pathways decisively different 
from conventional QD carrier dynamics, especially for holes, for which the piezoelectric side 
barriers (shown in Fig. [2]) block the relaxation to the DP minimum along the y axis. The 
relaxation of holes along the x axis to the DP minimum is hindered as well, but by the PEP 
minima which capture the holes before they can enter the DP minima. The inefficient hole 
relaxation was also observed in the experiments as weak QD luminescence during Ti-Sap 
carrier generation directly to the QW. This indicated a very inefficient carrier relaxation 
from the QW to the radiative states of the SIQD's unless the carriers of the QW were hot, 
as in the case of Ar-ion laser pumping. 
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On the basis of these observations we have developed the carrier dynamics model, 
schematically shown in Fig. [HI This model is a generalization, based on a true three- 
dimensional description of the QD system, of the models presented in Refs. 9| and [3], 
which are as such not able to describe the experiments analyzed here. 



A. Underlying assumptions 

Our carrier dynamics model is based on the following assumptions, which are motivated 
by both experimental observations and theoretical studies: 

(i) The carrier population is charge neutral. This implies that the total number of holes 
(including the valence states of the DP minimum, the two PEP minima and the QW) is 
equal to the total number of electrons (including the conduction states of the DP minimum 
and the QW). This is motivated for undoped samples were the optical carrier generation 
and recombination conserves the net charge, except for small charge fluctuations between 
different QD's, which mainly broaden the PL peaks. The assumption of charge neutrality 
is also supported by the invariance of the experimentally measured line energies in SIQD 
luminescence spectra, as a function of state filling; i.e., the luminescence energies are inde- 

. If there was a change in the net charge, there 
would also be a noticeable change in the line energies. 

(ii) The total number of electron and hole states, included in the model, were equal. The 
number of electron and hole states was limited to ~ 100, in order to stabilize the MC 
simulations and to account for the possible saturation of carrier generation to the QW 
reservoirs. This was done by setting the total number of electron states equal to the total 
number of hole states of the entire model. This has little effect on the QD luminescence 
character as we are typically far from saturation, but this is in accordance with assumption 

(i). 

(iii) We excluded the possibility of nonradiative recombination channels. Nonradiative 
recombinations are typically a result of crystal imperfections or interface states and they are 
indirectly observed in terms of very short luminescence lifetimes. The self-organized samp 
of SIQD's are of very high quality, showing very long luminescence life times of ~ Ins. 
The non-radiative recombination channels are thus rare. 

(iv) Nearly degenerate electron and hole states were combined to a few degenerate levels. 



es 
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The DOS had to be simphfied in order to make the MC calculations numerically feasible. 
The original DOS contained too many states to be included in the MC model. However, the 
simplified DOS resembles well the true multiband three-dimensional DOS. 

(v) We accounted for the spin of the carriers. The preservation of the spin quantum 
number infiuences the total speed of recombination. The conservation of the spin should 
be accounted for in any carrier dynamics model, however, it has to our knowledge not been 
accounted for in previous QD relaxation models. Note, that the spin is in general not a good 
quantum number for states of the eight-band effective-mass Hamiltonian, since the valence 
bands are strongly coupled by the spin-orbit interaction. In the present case, the spin is, 
however, an approximately good quantum number for both electron and hole states. This is 
due to the vertical confinement and biaxial strain of the QW, which cause the most strongly 
confined hole states to be nearly pure heavy-hole states, which have a well defined spin. 

(vi) The intrahand relaxation time Tq was assumed ten times smaller than the radiative 
life times. An equal intraband relaxation time of both electrons and holes is motivated for 
high carrier densities in which the carrier relaxation is dominated by Coulomb scattering and 
the Auger process. 3, 7| From experiments!^, Q] we know that at relatively high state filling 
the intraband relaxation is one order of magnitude faster than the radiative recombination. 
Furthermore, we assumed that relaxation takes place between consecutive states only. 

(vii) The tunneling rate between the PEP and DP minima was assumed zero. This tun- 
neling rate is small in comparison with the prevailing relaxation and recombination rates, 
since the spatial separation of the DP and PEP minima is about 20 ~ 30 nm. 

(viii) It was assumed that an intense THz radiation can couple the hole states pi and 
hi. If the intensity of the THz radiation is large and if the THz photon energy matches 
the energy separation of hi and pi it gives rise to a resonant charge transfer where the 
absorption (hole: pi — hi) and emission (hole: hi pi) of THz photons compete. The 
radiation field is assumed to be in resonance with the hi DP state and any of the deep 
PEP states, which are in the model described by a single degenerate level pi. We neglect 
spontaneous emission, which is presumably small at high THz intensities. The strength of 
the hole transfer is, in the experiment, not limited to a very narrow THz frequency range 
since in the real SIQD sample the single resonance frequency of the model is replaced by a 
quasi- continuum of transition frequencies. The inhomogeneous size distribution of the QD's 
is likely to further broaden the otherwise narrow peaks of allowed THz radiation frequencies 
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of effective coupling. 

The coupling strength between the PEP states pi and the DP states hi depends on 
the intensity of the THz field, the carrier populations of the involved hole levels and a 
phenomenological coupling constant. This coupling strength is the only free scaling constant 
in our model. The available experimental and numerical tools are currently far from enabling 
quantitative measurements or ab initio type calculations of this coupling constant. As a 
summarizing note on our relaxation model we recall that all model parameters, except the 
phenomenological THz radiation coupling, are based on direct measurements (r^j, tq), fitting 
to experimental data (/j), or electron structure calculations ((?«). j^, Q, Q 



B. Master equations 

The time evolution of the electron (raf ) and hole [n^) populations of the intermediate DP 
states are described by the following master equations for 2 < i < 4, 



dni <+i {Qi - n'r) n\ {gi^i - n^i) nfn'^ 



dt To To QiTr 



^>Um + ''l±})A^^_f^)^ (3) 



h (4) 



dni _ nil {g, - ) {g,^i - nl,) nln^ 

dt To Tq QiTri 

and for the PEP states with 2 < z < 5 

dt To To 

The time evolution of the ground state {i = 1) populations are given by 

dnl nl {gi - n^) n|n^ nl \n\j^ + n^g] 

"TT Jl 7. IDJ 

dt To giTri 2giTrl 

dn\ _ {gi - Q _ n^n^ ^ 

dt To grlTrl 

+ [gi{n{^ + nU-2n'yi]GTH., (7) 

drii^ ^27 {9l - ^17) , (^hv p p\ ^ /o^ 

— r [Uigi - rti^gi) (jcTHz, \o) 
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and for the reservoirs {i = R) we have 

dn^ ^ n% (c/4 - nl) _ n\n\ 

dt TeQD QrTr 



+ G, (9) 



dn),_ <(^74-<) n\{2gl-nl^-nl^) 



dt ThQD 'TpQD 



h 



- ^ + G, (10) 

where the two PEP minima are identified by subindices = iA and 7 = iB. In Eqs. IMTOl 
is the THz radiation intensity and G = {gfj^ — ra^) ((7^ — n^J Gexc is the carrier generation 
rate by the pump laser Gexc- We also accounted for the spin of the carriers, which effectively 
raised the population of excited states. A detailed analysis of Eqs. [BlfTUl follows below. 



C. Partial transition rates 



The MC carrier dynamics model is depicted in Fig. [H] and was given in mathematical form 
in Eqs. l3llT0l It describes the time evolution of the electron and hole populations of 
the DP states as well as the hole populations nf^ and n^^ of the PEP states (left and right 
minima) . 



1. The population of intermediate states 

The time evolution of the electron and hole populations and of the intermediate 
states (i.e., 2 < i < 4) of the DP minima are described by Eqs. [HJand |H The population of 
state i is increased by the intraband relaxation (colored/gray arrows in Fig. ^ from state 
z + 1 to state i [first term on the right hand side (RHS) of Eqs. [3]and|l] and reduced by the 
carrier relaxation to state i — 1 [second term on the RHS of Eqs. [3]and|4]. The rate of the 
intraband relaxation is described by the intraband relaxation lifetime tq, which is typically 
one tenth of the radiative lifetime r^.jsl 

The main recombination processes involve electron and holes (the center black dashed 
arrow in Fig. [H]), both confined in the DP minima, and are described by the third term on 
the RHS of Eqs. [3] and HI The recombination of electrons in the DP minimum and holes 
in the PEP minima is described by the fourth term on the RHS of Eq. [3] (left and right, 
black dashed arrows in Fig. [6]). This rate is weaker than the recombination rate between the 
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DP states, due to the small spatial overlap of the pertinent wave functions. The different 
transition probabilities of the DP^DP and DP^PEP recombination could be included 
by using different time constants. However, in our model it has been more convenient to 
describe this effect by the ratios — fi)- The numerical values of fi were determined 
by fitting the state fil ling (P L peak intensities), as a function of the generation rate, to 
experimental PL data. 12l. Il3| 

The time evolution of the hole populations of the intermediate (i.e. 2 < z < 5) states of 
the PEP minima are described by Eq. O The first and second terms on the RHS correspond 
to intraband relaxation (colored/gray arrows in Fig. [6]) from state i + 1 to state i and from 
state z to z — 1, respectively. The last term corresponds to radiative recombination with 
electrons confined in the DP minimum (left and right, black dashed arrows in Fig. [6]). 



2. The population of the DP and PEP ground states 

The time evolution of the electron ground state population nl is described by Eq. [6], which 
is similar to Eq. [HI except that it does not contain the second term of the RHS of Eq. [H] as 
the holes cannot relaxate any further from the ground state. Equation [7] describes the time 
evolution of the hole population of the DP ground state. The first and second terms on 
the RHS correspond to intraband relaxation from state i + 1 to state i (colored/gray arrows 
in Fig. EI) and to radiative recombination (center black dashed arrow in Fig. [6]) between 
electron and holes of the DP minima. The third term correspond to the THz radiation- 
induced coupling of the states hi, plA (in the left PEP minimum) and plB (in the right 
PEP minimum). It accounts for the absorption of a THz photon, exciting a hole from state 
hi to states plA or plB, and the stimulated emission of a THz photon, taking a hole from 
state pi A to hi or from state plB to hi. These processes are indicated by double-headed 
arrows in Fig. El 

The time evolution of population nf^ of the deepest PEP states is described by Eq. |8l 
where 7 = A or 7 = i? correspond to left and right PEP minimum, respectively. The first 
term on the RHS corresponds to the intraband relaxation from state p2 to state pi (col- 
ored/gray arrows in Fig. [6]). The second term accounts for both the emission or absorbtion 
of a THz photon (simultaneously moving a hole from state pl-f to state hi). 
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3. The population of the reservoirs 

The physical QW, embedding the SIQD's is modelled as a reservoir of electrons and holes, 
with the state degeneracies (^Ij and both being much larger than the degeneracies of the 
QD levels. The electron and hole reservoirs are described by Eqs. [9] and [101 respectively. 
The first term on the RHS of Eq. [9] [Eq. [10] describes the capturing of an electron (hole) 
from the reservoirs to the DP minimum with the capture lifetime TeQ/j (t/iqz?) describing the 
rate of this process. The second term on the RHS of Eq. [HI and the third term on the RHS of 
Eq. [To] account for the radiative recombination of electrons and holes in the reservoirs, with 
the corresponding lifetime tr. The second term of Eq. [10] is due to holes getting captured 
into the two {A and B) PEP minima. The carrier generation is described by the last term 
of both Eq. [9] and [10] 



D. General remarks on the model parametrization 

We note that the radiative lifetimes of SIQD's have been measured and modeled suc- 
cessfully previously. In our carrier dynamics model we used radiative lifetimes T^i 
{Tri > Trj ~ 10 Tq, for all i > j), which were calculated from the eight-band electron 
structure model in the electric dipole approximation and were in good agreement with ex- 
periments. The radiative recombination rates of the pi states were found negligible. We 
found that the radiative recombination rates of also the lower excited PEP states {pi, with 
i < 4) were noticeable longer than those of the corresponding DP states. We used, there- 
fore, recombination lifetimes = Tri/{1 — fi) with < /j < 1 for the excited states of the 
PEP minima. The exact values of were determined by fitting the state filling (PL peak 



intensities) to experimental data. 13|] The distinction between DP and PEP states becomes, 
however, difficult for highly excited states, since these states are not confined in only one 
minimum, but the probability densities are shared between all three minima. 

The intraband relaxation time constant tq affects, together with the amplitudes of the 
individual recombination processes, the relative CW intensities and the rates of fading off 
of the PL peaks. The effect of the THz radiation on the QD PL is not very sensitive to tq 
at low generation intensities. However, the smaller tq is, the smaller generation intensity is 
needed to saturate the hole population of hi (even without any THz radiation). The value 
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of To does thereby affect the threshold value of the generation intensity at which the THz 
radiation-induced PL effect disappears. 

The validity of the used radiative and intraband relaxation time constants, describing the 
carrier dynamics in the absence of any THz radiation, was verified by comparing the numer- 
ical time dependent PL of our model with experiments. On the basis of this comparison we 
conclude that our model is very well in line with Refs. 13| and 7 

Equations [BlfTUl were solved by time-dependent MC simulations 
independent QD's. We obtained the probabilities of each many-particle partition as a func- 
tion of time by averaging over the ensemble populations. By simulating QD ensembles large 
enough (a few thousands of QD's) we were able to minimize the time-dependent oscillations 
of simulated quasistationary partitions during steady state laser modulation. 



(see also Fig. E]). 



261 ] of a large ensemble of 



IV. DISCUSSION AND CONCLUSIONS 

A. Luminescence during continuous pumping 

The free carriers can relax from the GaAs barrier to the DP and PEP minima, either 
directly or through the QW state continuum. The direct relaxation is assisted by the funnel 



shaped deformation and piezoelectric potential minima in the barriers. 22|] As the holes are 
predominantly confined in the deep PEP minima, the radiative recombination starts only 
when these minima are already filled with holes. A saturated QD is consequently strongly 
polarized during continuous generation and recombination. The charge separation, with 
holes in the PEP minima and electrons in the DP minimum, can persist even seconds after 
turning off the carrier generation. 

The picture described above does not account for any excitonic carrier-carrier interactions. 
Previous studies have, however, shown that the carrier-carrier correlation effects are small in 
SIQD's. 4] The very good qualitative agreement between the experimental CW PL and the 
PL simulated here does also support the used single-particle approximation. The exciton 
effects of SIQD's are, though, fairly vaguely understood and do call for further studies. 
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B. Influence of THz radiation on the steady-state luminescence 

According to our simulations, the enhancement of the QDl luminescence is due to a 
THz radiation induced continuous drift of holes between the PEP pi levels and the DP hi 
level (Fig. [8]). Under steady-state condition the total integrated PL intensity is conserved 
in the first approximation. The QD2, QD3, QD4, and QW peaks are accordingly reduced 
when the QDl peak is increased by the THz radiation, both in the experiment and theory. 
Our calculations predict furthermore that the THz radiation enhancement of the QDl peak 
disappears at high Ar+ pumping intensities as a result of the saturation of the hole population 
at hi. 

The results are in good qualitative agreement with the experiments, although, a quan- 
titative agreement cannot be achieved with the current model. Our model does not show 
as dramatic quenching of the QW PL by the THz radiation as was seen in the experi- 
ment, although. Figs. [7|and[H]do show a clear decrease of the QW PL. We argue that the 
large decrease of the experimental QW PL is related to the heating and ionization of QW 
excitons. [sl 

The discrepancies between the intensity distributions of the PL peaks (in Fig. [7j) of 
our model and of the experiments are mainly due to the simplified electronic structure 
used in our MC model and the fact that we have used the same relaxation life time for 
all hole and electron states. A better fit between the simulations and experiments could 
have been obtained by using, e.g., several different relaxation rates between the different 
hole levels. This would, however, have made the model less transparent. The exact PL 
intensity distribution between the different peaks is also very much dependent on the carrier 
generation rate. This is true for both the simulations and the experiments (see, e.g., Ref. j^). 
Figure 4 was plotted for a carrier generation rate where the simulated QD2 peak already 
exceeds the QDl peak. This generation rate was chosen, in order to obtain as clear THz 
effect as possible. 

C. THz radiation-induced transients 

Figures [H](a) and[8](b) show experimental and simulated time-resolved PL. The black and 
red curves correspond to integrated PL of the QDl and QW peaks, respectively. The left 
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panels show results obtained for overlapping time windows of Ti-Sap pumping and THz 
radiation excitation. The QW peak is decreased and the QDl PL is increased during the 
THz radiation with both PL peaks returning to their initial value after turning off the THz 
radiation. Note the saturation of both the experimental and theoretical PL during THz 
radiation. 



D. THz radiation-induced delayed PL 

The right panels of Fig. [8] show PL with a delayed THz radiation pulse after turning 
off the pump laser. The rise of the QDl peak by the delayed THz radiation is due to the 
excitation and drift of trapped holes from the PEP minima to the DP minimum, where 
they recombine with electrons localized in the DP minimum. The simulation time was much 
smaller than that of the experimental time scale, due to computational limitations. However, 
the time scale of the simulations were chosen long enough to observe saturation of the PL 
both during the pumping and between the pumping and the delayed THz radiation. There 
is also in this case a very good qualitative agreement between experiments and theory. 



E. Conclusions 



We have shown that the anomalous cooling of radiative QD excitons[29l] is due to a THz 
radiation-induced drift of holes from piezo-electric potential (PEP) minima to the lowest 
radiative states of the QD. The subsystem of radiative excitons is cooled, although, the 
average energy of the total carrier system increases (heated). 

By cooling of the radiative exciton subsystem we mean that the average hole population 
of the ground state hi increases whereas the average population of all excited states in the 
deformation potential (DP) minimum decreases. The large shift of the particle population 
(within the radiative states) from the excited states towards the ground state is analogous to 
lowering of the temperature in a system of thermal equilibrium. However, the carriers of our 
SIQD's are not in thermal equilibrium nor in thermal quasi-equilibrium during continuous (or 
transient) carrier generation and THz modulation. Therefore, the reduction of the average 
energy of the radiative excitons cannot be called lowering of the carrier temperature. 

Our computational model reproduces currently the following experimental observations 
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on SIQD's. (i) The simulated CW PL reproduces experimental data[l3| as a function of the 
carrier generation intensity, (ii) The THz radiation increases the ground state luminescence 
and decreases the excited state luminescence at low generation intensities during CW THz 
excitation, (iii) At the onset of the THz radiation there is a strong and sudden increase of the 
QDl luminescence. The enhancement by the THz radiation is then, however, exponentially 
reduced and saturates at a lower level, (iv) The THz radiation gives rise to ground state 
luminescence long after turning off the carrier generation. This delayed PL flash is seen only 
in the ground state luminescence and it fades out exponentially. 

The ultimate goal of our work is to motivate new experiments in which the THz photon 
energy would be varied over a wide energy range, covering several pi and hi levels. We 
argue that the analysis of the original experiments by Yusa et al. is necessary for a careful 
planning of further experiments. Finally, we note that our simulations suggest that the 
resonance transfer by THz radiation could be replaced by a bias voltage-driven resonance 
tunneling, leading to an enhanced QD ground state luminescence. 
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TABLE I: Material parameters from Landolt-Bornstein (Ref. [j]) unless otherwise noted 



Parameter 


GaAs In 


o.iGao.gAs 


InP 


flo (A) 


5.65 


5.69 


5.87 


cii (10^-*^ dyn cm~^) 


11.90 


11.50 


10.11 


Ci2 (10^-*^ dyn cm^^) 


5.38 


5.29 


5.61 


C44 (lO^-"^ dyn cm~^) 


5.95 


5.75 


4.56 


ei4 (C m-2) 


-0.16 


-0.15 


-0.04 (Ref. [19]) 




12.53 


12.70 


13.90 



TABLE IL Eight-band k • p parameters from Ref. 25 1 



Parameter 


GaAs 


Ino.1Gao.9As 


Eg (eV) 


1.519 


1.366 


A (eV) 


0.341 


0.332 


71 


2.556 


3.464 


72 


-0.152 


0.295 


73 


0.718 


1.148 




0.067 


0.062 


Ep'' (eV) 


20.160 


19.742 


VBO (eV) 


-0.055 


0.000 


a„ (eV) 


1.160 


1.144 


ac (eV) 


-7.170 


-7.196 


by (eV) 


-2.000 


-1.980 


dy (eV) 


-4.800 


-4.680 




12.530 


12.727 



"These are the used value of Ep and are slightly smaller than the values given in Ref. [25|. This was as 



explain^ in the text. 
^Ref. 
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(a) QD Sample 

-Ar+ or Ti-Sap » 





291]) the QD sample was pumped 



FIG. 1: (Color online) (a) In the analyzed experiment (Ref. 
by either an Ar+ or Ti-Sap laser. The QD intraband dynamics was modulated by a free electron 
laser (FEL). (b) Strain-induced quantum dot in an Ino.2Gao.8As QW, strained by an InP island. 



The island is defined by the crystal planes {001}, {101}, and {111}. 
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FIG. 2: (Color online) Potential energy of (a) electrons and (b) holes, in the middle of the QW 
(8.25 nm below the InP island). The energies are given with respect to the band edges of an 
Ino.2Gao,8As QW, far away from the island. The white dashed lines show the bottom contour of 
the InP island and the labels correspond to the local potential maxima/minima. 
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FIG. 3: (Color online) (a) Electron and (b) hole probability densities of the four eigenstates, lowest 
in energy. The position of the QW is shown by the gray faces. 
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FIG. 4: (Color online) Density of hole and electron states, including states confined in the DP and 
PEP minima. The blue/gray curve is the total DOS while the black curve corresponds to the DOS 
of carriers confined in the DP minima only, i.e., the optically active states. The top panels (a)-(f) 
show the probability densities of selected states, with respect to the position of the InP island (red 
lines) above the QD. 
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FIG. 5:^ (Color online) Simulated CW PL (gray/colored lines) in comparison with experimental 

(Ref. [1^) data (black solid line). Figure [5](a) shows the effect of the InP island geometry. 
Both theoretical PL curves are based on full 3D simulations, including the piezoelectric potential. 
However, the solid blue/gray line was computed with InP islands having the shape of a truncated 
cone (axisymmetric), whereas the dashed line correspond to simulations based on the InP island 
geometry of Fig. [H Figure [5](b) shows the effect of the piezoelectric potential in the case of an 
axisymmetric InP island. The numerical PL curves have been normalized by equating their QD4 
peak to the corresponding QD4 peak of the experiments 
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FIG. 6: (Color online) Carrier dynamics model. The eigenstates of the DP minima are labeled 
Cj and hi] the hole states of the PEP minima are labeled pi. The pump laser G generates carriers 
to the reservoirs {cr and /i/j) and the resonant THz radiation transfers holes between the ground 
states pi and hi . Radiative recombinations (black dashed arrows) are shown for the third electron 
level only. The energy scale is only schematic. 
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FIG. 7: Experimental (a) and theoretical (b) photoluminescence from a sample pumped by an 
Ar''" laser only (solid lines) and by Ar"*" laser and FEL simultaneously (dashed line). The FEL 
energy of the experiments was Hiothz = 10.4 meV. 
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FIG. 8: (Color online) Experimental (a) and simulated (b) time-resolved PL of the QD ground 
state (black) and QW (red) for T = 15 K. The energy and power of the THz radiation were 
Hlothz = 2.5 meV and P = 1 kW, respectively. The time windows of the Ti-Sap laser (carrier 
generation to the QW) and THz radiation are indicated by horizontal bars. 
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